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Abstract 

Background: Downy mildew, caused by Plosmoporo viticolo, is one of the most severe diseases of grapevine and is 
commonly controlled by fungicide treatments. The beneficial microorganism Trichoderma horzionum T39 (T39) can 
induce resistance to downy mildew, although the molecular events associated with this process have not yet been 
elucidated in grapevine. A next generation RNA sequencing (RNA-Seq) approach was used to study global 
transcriptional changes associated with resistance induced by T39 in Vitis vinifera Pinot Noir leaves. The long-term 
aim was to develop strategies to optimize the use of this agent for downy mildew control. 

Results: More than 14.8 million paired-end reads were obtained for each biological replicate of T39-treated and 
control leaf samples collected before and 24 h after P. viticola inoculation. RNA-Seq analysis resulted in the 
identification of 7,024 differentially expressed genes, highlighting the complex transcriptional reprogramming of 
grapevine leaves during resistance induction and in response to pathogen inoculation. Our data show that T39 has 
a dual effect: it directly modulates genes related to the microbial recognition machinery, and it enhances the 
expression of defence-related processes after pathogen inoculation. Whereas several genes were commonly 
affected by P. viticola in control and T39-treated plants, opposing modulation of genes related to responses to 
stress and protein metabolism was found. T39-induced resistance partially inhibited some disease-related processes 
and specifically activated defence responses after P. viticola inoculation, causing a significant reduction of downy 
mildew symptoms. 

Conclusions: The global transcriptional analysis revealed that defence processes known to be implicated in the 
reaction of resistant genotypes to downy mildew were partially activated by T39-induced resistance in susceptible 
grapevines. Genes identified in this work are an important source of markers for selecting novel resistance inducers 
and for the analysis of environmental conditions that might affect induced resistance mechanisms. 
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Background 

Plasmopara viticola (Berk, and Curt.) Berl. and de Toni 
is a biotrophic oomycete that causes downy mildew in 
grapevine [1]. This devastating disease occurs worldwide, 
particularly in regions with warm and wet conditions 
during the growing season. P, viticola mainly infects 
leaves and clusters of young berries and produces oil 
spot lesions on the adaxial leaf surface accompanied by 
massive sporulation on the abaxial surface. Although 
downy mildew can be controlled by frequent applications 
of chemical fungicides, concerns about the environmen- 
tal impact of pesticide overuse [2] and the development 
of resistant P. viticola populations [3] have sparked an 
interest in alternatives to chemical treatments. 

The grapevine industry relies predominantly on Vitis 
vinifera, which is susceptible to downy mildew. Resist- 
ance traits have been identified in wild grapevine spe- 
cies, and the mechanisms of resistance to downy mildew 
have been characterized in resistant genotypes [1]. Tran- 
scriptomic analysis supports the view that downy mildew 
resistance is mainly a post-infection phenomenon [4] 
and highlights the importance of transcriptional repro- 
gramming in both resistant and susceptible genotypes in 
response to P. viticola inoculation [4-8]. Transcriptional 
changes associated with P, viticola infection of suscep- 
tible grapevines have been related to a weak defence re- 
sponse [4] and to the establishment of a compatible 
interaction [5,7,9,10]. The response of resistant geno- 
types has been characterized by strong and rapid tran- 
scriptional reprogramming of processes related to 
defence, signal transduction, and secondary metabolism, 
which are either not induced or induced to a lesser ex- 
tent in susceptible grapevines [4,8,11-14]. In particular, 
downy mildew resistance has been correlated with 
enhanced expression of genes encoding pathogenesis- 
related (PR) proteins and enzymes of phenylpropanoid 
biosynthesis, and with specific modulation of signal 
transduction components and markers of hypersensitive 
response (HR) in resistant grapevines [4,8,11-13]. 

Downy mildew symptoms can be significantly reduced 
in susceptible grapevines by preventive application of 
resistance inducers, such as chitosan [15], laminarins 
[16-18], p-aminobutyric acid (BABA) [19,20], acibenzolar- 
S-methyl (BTH) [21,22], and thiamine [23]. Treatments 
with plant extracts [24] or microbial extracts [22,25] have 
also been found to increase grapevine resistance to downy 
mildew. The ascomycete Trichoderma harzianum strain 
T39 (T39) significantly reduces downy mildew symptoms 
by activating grapevine resistance both locally and system- 
ically [21], although the molecular events responsible for 
resistance induction have not yet been clarified. Tricho- 
derma spp. strains have been characterized in model sys- 
tems based on their ability to induce plant resistance 
against pathogens [26,27] by reprogramming the plant 



transcriptome [28-32]. Specific strains of beneficial micro- 
organisms can improve plant performance by activating a 
plant-mediated defence mechanism known as induced 
systemic resistance (ISR) [33]. Through root or leaf inter- 
actions [21,34], beneficial microorganisms are recognized 
by the plant, which results in a mild but effective activa- 
tion of the plant immune responses in all tissues [35]. ISR 
confers broad-spectrum resistance to various types of 
pathogens and abiotic stresses [27,36] and is usually regu- 
lated by jasmonic acid (JA)- and ethylene (ET) -dependent 
signalling pathways [33]. Rather than directly activating 
plant defences, beneficial microorganisms that induce re- 
sistance usually prime the plant so that it responds more 
strongly upon exposure to the stress condition [35,36]. 
Primed plants display faster and/or stronger activation of 
the defence responses after pathogen inoculation [37]. Be- 
cause plant defences are activated only when they are 
really needed [38], priming provides advantages in terms 
of energy costs for the plant [39], and it is probably 
evolved to save energy under pathogen-free conditions 
[40]. The benefits of priming outweigh its costs when dis- 
ease does occur, and priming is seen as a promising strat- 
egy in modern disease management [41]. The absence of 
apparent energy costs associated with T39-induced resist- 
ance in grapevine suggested a priming state activation 
[42]. However, the molecular mechanisms underlying re- 
sistance induction have been only partially revealed 
[26,43], and information regarding induced resistance in 
grapevine is particularly scarce [20,44]. 

In this study, we used Illumina RNA-Seq analysis to 
characterize the global transcriptional dynamics associated 
with T39-induced resistance in grapevine. To the best of 
our knowledge, this study is the first to use high- 
throughput sequencing technology to investigate molecu- 
lar events underlying induced resistance in plants, and it is 
also the first transcriptome-wide characterization of resist- 
ance induced by a beneficial microbe in a non-model 
plant. Our analysis revealed that the reduction of downy 
mildew symptoms is related to a complex transcriptional 
reprogramming in T39-treated plants, both before and 
24 h after pathogen inoculation. In particular, the reaction 
of T39-treated plants to pathogen inoculation is associated 
with enhanced expression of P, viticola-responsive genes 
and specific modulation of some genes related to defence 
in resistant grapevines. Our study has identified genes that 
could be valuable as markers of ISR activation for subse- 
quent selection of new resistance inducers with improved 
ability to stimulate plant defences. 

Results 

RNA-Seq sequencing and mapping of reads to the 
grapevine genome 

The ability of T39 to activate local and systemic resist- 
ance processes has been previously reported in grapevine 
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[21,42], and leaf treatment was used to study responses 
to T39 treatment and pathogen inoculation in the same 
plant organ. Leaf samples were collected before and 24 h 
after P, viticola inoculation from T39-treated and con- 
trol plants, and resistance induction was confirmed by 
phenotypic analysis (Figure 1). Four different treatments 
were analysed by RNA-Seq: control (C), T39-treated 
(T39), viticolaAnoculdited control (C+P.v.) and P, viti- 
co/<3^-inoculated T39-treated (T39+P.V.) plants. Three 
biological replicates of each treatment were analysed, 
with each biological replicate sequenced twice in separ- 
ate lanes (sequencing replicate). Between 6.1 and 18.8 
million paired-end reads of 100 nucleotides were 
obtained for each sequencing replicate, and an average 
of 94% of these passed the quality control test (filtered 
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Figure 1 Grapevine resistance induced by treatment with 
Trichoderma harzianum T39. Downy mildew severity (%) was 
assessed on Vitis vinifero Pinot Noir control plants (C) and plants 
treated with T. harzianum T39 (T39). Mean severity scores and 
standard errors for six replicates are presented. Treatments followed 
by different letters are significantly different according to Tukey's 
test {P < 0.05). Photographs show sporulation of downy mildew on 
representative leaves. 

V J 



reads; Additional files 1 and 2). Summing the reads of 
two sequencing replicates, more than 14.8 million fil- 
tered reads were obtained for each biological replicate 
(Table 1), corresponding to a coverage of at least 32 x 
the grapevine transcriptome (Table 2). 

An average of 91% of filtered reads mapped to the 
grapevine genome (mapped reads) for each biological 
replicate (Table 1), and similar percentages mapped for 
each sequencing replicate (Additional file 1). An average 
of 96% of mapped reads matched to unique locations 
(unique reads), 4% displayed multiple matches (multi- 
reads) to the grapevine genome (Table 1 and Additional 
files 1 and 3), and about 77% of mapped reads matched 
to grapevine genes (Additional files 3 and 4). 

Gene expression estimation by RNA-Seq 

Gene expression level was assessed on the basis of both 
unique and multi-reads (Additional file 3) to improve 
evaluation of members of multigene families [45,46], and 
it was expressed as fragments per kilobase of transcript 
per million fragments mapped (FPKM) values using Cuf- 
flinks [46]. Between one and nine transcript isoforms 
were recognized for each grapevine gene, and 3,679 to 
7,548 novel genes were identified in each sequencing 
replicate (Additional file 5), representing important new 
information for genome annotation. About 66% of 
grapevine genes were expressed in each sequencing rep- 
licate (Additional files 5 and 6). High correlations be- 
tween sequencing replicates were obtained (Pearsons 
correlation coefficients and values greater than 0.98 
and 0.95, respectively; Additional files 7, 8 and 9), and 
read counts were summed to obtain better coverage and 
improve variance estimation of each biological replicate 
[47]. Considering the expression of all grapevine genes, 
the Pearson s correlation coefficient between control and 
T39-treated plants (T39 vs, C) was greater than that be- 
tween control and P, vitico la-inoculated plants (C+P.v. 
vs, C and T39+P.V. vs, C; Table 3), suggesting that few 
transcriptional changes were caused by T39 treatment 
and that major transcriptional reprogramming occurred 
after pathogen inoculation. 

Grapevine genes differentially expressed during 
Trichoderma harzianum T39-induced resistance 

DeSeq statistical analysis [47] revealed that 7,024 genes 
were differentially expressed with respect to the control 
with a false discovery rate (FDR) of 5% and a minimum 
fold-change of two in at least one pairwise comparison 
(Additional file 10). About 90% of these genes were 
expressed in all treatments (Figure 2a), indicating that 
T39 and P, viticola caused a significant modulation (up- 
or down-regulation greater than 2-fold) of genes nor- 
mally expressed in control leaves. Relatively few genes 
were specifically expressed in response to the T39 
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Table 1 RNA-Seq sequencing and read mapping to the grapevine genome 



Treatment^ 


Replicate'' 


Total reads'^ 


Filtered reads'^ 


% 


Mapped reads^ 


% 


II- J f 
Unique reads 


% 


Multi-reads^ 


% 


C 


1 


20480361 


1 9388264 


95 


1 7975427 


93 


17473619 


97 


501808 


3 




2 


20182309 


19095859 


95 


17935380 


94 


1 7392482 


97 


542898 


3 




3 


15856292 


14845528 


94 


13327853 


90 


12749078 


96 


578775 


4 


T39 


1 


24312739 


23208015 


95 


21 103884 


91 


20207496 


96 


896388 


4 




z 


DDD \ j4Uo 


OZJOJOZD 


yz 


oU4UyUoD 




zy4yjo 1 z 


y/ 


91 3274 


T 

D 




3 


25081653 


23722385 


95 


21873968 


92 


21145558 


97 


728410 


3 


C+P.v. 


1 


21934228 


20800485 


95 


19240139 


92 


18603657 


97 


636482 


3 




2 


25097152 


23684317 


94 


21637677 


91 


20896985 


97 


740692 


3 




3 


22471334 


21555054 


96 


19384445 


90 


18634940 


96 


749505 


4 


T39+P.\/. 


1 


22897386 


21318685 


93 


18406407 


86 


17576560 


95 


829847 


5 




2 


29049603 


26710912 


92 


24801196 


93 


23981057 


97 


820139 


3 




3 


22603588 


21369934 


95 


19558774 


92 


18859902 


96 


698872 


4 



^ Grapevine leaves of control (C), Trichoderma harzianum T39-treated (T39), Plasmopara wf/co/a-inoculated control {C+P.v.) and P. Wf/co/a-inoculated T39-treated 
(T39+P.V.) plants. 

^ Biological replicates (plants), numbered from 1 to 3. 

Total sequenced reads as sum of reads obtained from two replicates of sequencing for each treatment. 
^ Reads passing the quality check and the corresponding percentage (%) of total reads. 

^ Reads mapping to Pinot Noir grapevine genome [77] Release 3 [78], and the corresponding percentage (%) of filtered reads. 
^ Reads mapping to unique locations in the grapevine genome, and the corresponding percentage (%) of mapped reads. 

^ Reads mapping to more than one location (2-100 matches) in the grapevine genome, and the corresponding percentage (%) of mapped reads. 



treatment (12 genes) and after P, viticola inoculation in 
control (36 genes) and T39-treated plants (22 genes). 

The expression profiles of the differentially expressed 
genes were determined by cluster analysis, and genes 
were grouped into 10 clusters based on their expression 
modulation (Figure 2b and Additional file 10). Among 



Table 2 RNA-Seq sequencing and coverage of the 
grapevine transcriptome for each biological replicate 



Treatment^ 


Replicate^ 


Sequenced 
bases (Mbp)" 


Read length 
(bp)^ 


Coverage 
(fold)^ 


C 


1 


1844 


95 


42 




2 


1812 


95 


42 




3 


1393 


93 


32 


T39 


1 


2195 


94 


50 




2 


3054 


94 


70 




3 


2253 


95 


52 


C+P.v. 


1 


1981 


95 


46 




2 


2244 


95 


52 




3 


2049 


95 


47 


T39+R\/. 


1 


1999 


94 


46 




2 


2500 


94 


57 




3 


2028 


95 


47 



^ Grapevine leaves of control (C), Trichoderma harzianum T39-treated (T39), 
Plasmopara wf/co/a-inoculated control [C+P.v) and P. v/f/co/a-inoculated 
T39-treated (T39+P.V.) plants. 

^ Biological replicates (plants), numbered from 1 to 3. 

^ Total bases (Mbp) sequenced by RNA-Seq analysis that passed the quality 

check. 

Mean length (bp) of sequenced reads. 
^ Coverage of the Pinot Noir grapevine transcriptome [77] Release 
3 (43.5 Mbp) [78]. 



genes modulated by T39, 343 had comparable expres- 
sion levels before and after P, viticola inoculation in 
T39-treated plants (cluster 1), while 233 genes were ini- 
tially modulated by T39 and then further modulated by 
P, viticola inoculation in T39-treated plants (cluster 2), 
In addition, the expression of 479 genes was affected ex- 
clusively by T39 treatment and not by P, viticola (cluster 
3). P, viticola inoculation resulted in modulation of 
3,454 genes with comparable expression levels in control 
and T39-treated plants at 24 h after inoculation; of 
these, 948 genes were directly modulated by T39 treat- 
ment (cluster 4) and 2,506 genes were not (cluster 5). 
Interestingly, 868 genes showed reinforced modulation 
in T39-treated plants compared to control plants after P, 
viticola inoculation (ISR-primed genes), indicating 



Table 3 Pearson's correlation coefficients between 
grapevine treatments 



Comparison 


Pearson's correlation 




coefficient 


C vs. T39 


0.93 


C vs. C+P.v. 


0.83 


C vs. T39+P.I/. 


0.80 


T39 vs. C+P.v. 


0.80 


T39 vs. J^9+P.v. 


0.81 


C+P.v. vs. J?>9+P.v. 


0.92 



Pearson's correlation coefficients are based on the expression values of all 
grapevine genes, expressed as fragments per kilobase of transcript per million 
fragments mapped (FPKM), in control (C), Trichoderma harzianum T39-treated 
(T39), Plasmopara wf/co/c7- inoculated control [C+P.v), and P. Wf/co/a-inoculated 
T39-treated (T39+P.V.) plants. 
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enhancement of the grapevine defence reaction to the 
pathogen. Among these genes, 518 and 350 genes were 
modulated (cluster 6) or not (cluster 7) by T39 treat- 
ment, respectively. A total of 888 genes were modulated 
by P, viticola exclusively in control plants; of these, 267 
genes were directly modulated by T39 treatment (cluster 
8) and 621 genes were not (cluster 9). Conversely, 759 
genes were modulated by P, viticola exclusively in T39- 
treated plants (cluster 10, ISR-responsive specific genes), 
and they represent the specific reaction to P, viticola 
of T39-treated plants. ISR-responsive specific genes 



Figure 2 Grapevine genes differentially expressed during 
Trichoderma harzianum T39-induced resistance to Plasmopara 
viticola. (A) Venn diagram summarising the distribution of 7,024 
differentially expressed genes in control (C), T39-treated (T39), 
Plasmopara wt/co/c/-inoculated control {C+P.v) and P. viticola- 
inoculated T39-treated (J39+P.v) plants. Differentially expressed 
genes were identified by the DESeq package [47] with a false 
discovery rate (FDR) of 5% and a fold-change greater than two in at 
least one pairwise comparison. (B) Grouping of differentially 
expressed genes into 10 clusters based on the expression profiles. 
Genes modulated by T39 were distinguished as those with 
comparable expression levels before and after P. viticola inoculation 
in T39-treated plants (cluster 1), genes initially modulated by T39 
and further modulated by P. viticola inoculation in T39-treated plants 
(cluster 2), and genes affected exclusively by T39 (cluster 3). Genes 
modulated by P. viticola with comparable expression levels in 
control and T39-treated plants were classified as those with (cluster 
4) and without (cluster 5) direct modulation by T39. Genes 
modulated by P. viticola with greater modulation (> 1.5-fold) in 
T39-treated than in control plants (ISR-primed genes) where 
classified as those with (cluster 6) and without (cluster 7) modulation 
by T39 treatment. Genes modulated by P. viticola exclusively in 
control plants were classified as those with (cluster 8) and without 
(cluster 9) modulation by T39, whereas genes modulated by P. 
viticola exclusively in T39-treated plants were grouped in cluster 10 
(ISR-responsive specific genes). 



were mainly induced (63%), while genes modulated by 
P. viticola exclusively in control plants were mainly 
repressed (55%). 

Functional annotation of differentially expressed genes 

Differentially expressed genes were automatically anno- 
tated and then grouped into 15 functional categories of 
Gene Ontology (GO) biological processes (Figure 3); 
genes that could not be associated with any biological 
process category were assigned to the GO root (bio- 
logical process). Although the group of genes with un- 
known function was the largest (35%), 11 categories 
were represented to different degrees in differentially 
expressed genes compared to the grapevine transcrip- 
tome. DNA (5.0%) and protein (8.4%) metabolic pro- 
cesses were significantly underrepresented among the 
differentially expressed genes, but response to stress 
(5.9%), response to stimulus (9.8%), and signal transduc- 
tion (2.7%) categories were significantly overrepresented. 
Likewise, transport (6.2%) and metabolism of carbohy- 
drates (3.4%), lipids (2.9%), secondary compounds 
(1.1%), and energy (1.2%) were significantly overrepre- 
sented. The large fraction of differentially expressed 
genes with unknown function highlighted the identifica- 
tion of novel processes potentially relevant for the in- 
duction of plant resistance and for the response of 
grapevine to P. viticola inoculation. However, additional 
studies are required to better characterize the relevance 
of these genes in the mechanisms of T39-indiuced 
resistance. 
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Figure 3 Annotation of differentially expressed genes. Distribution of differentially expressed genes (black) and grapevine transcriptome 
(grey) among the 15 selected Gene Ontology (GO) functional categories (GO identifiers are reported in brackets). Frequencies are calculated as 
the percentage of the total number of GO biological process terms (8,407 for differentially expressed genes and 39,771 for the grapevine 
transcriptome) obtained using the ARG0T2 function prediction tool [83,84]. Asterisks indicate GO functional categories differentially represented 
in the differentially expressed genes compared to the entire grapevine transcriptome according to GOstat statistical analysis (P< 0.05) [87]. 



After genes with unknown function were discounted, a 
specific distribution of up- and down-regulated func- 
tional categories was observed in the various clusters 
(Figure 4). Genes of cluster 1 involved in carbohydrate, 
protein, and secondary metabolism were mainly induced 
by T39 (69, 76 and 83%, respectively), whereas those 
related to lipid metabolism and transport were mainly 
repressed (86 and 65%, respectively; Figure 4a). Genes of 
cluster 2 related to protein metabolism and transcription 
were mainly repressed by T39 and subsequently induced 
by P, viticola inoculation of T39-treated plants. More 
specifically, T39 treatment directly induced the expres- 
sion of genes encoding protein kinases and receptor-like 
protein kinases, MYB transcription factors, hormone- 
responsive genes, defence-related proteins, and enzymes 
of secondary metabolism. Annotation of genes modu- 
lated exclusively by T39 (cluster 3) revealed down- 
regulation of processes related to response to stimulus, 
cellular component organization, and protein metabol- 
ism (Figure 4b). 

Genes that were similarly modulated by P. viticola in 
control and T39-treated plants (clusters 4 and 5) were 
mainly related to primary metabolic processes and signal 
transduction (Figure 4c). In particular, genes involved in 
cellular component organization and carbohydrate me- 
tabolism were mainly induced (74 and 65%, respectively). 
More specifically, P. viticola induced six cellulose 
synthase-like genes, six expansin genes, as well as genes 
encoding enzymes of glucan biosynthesis and hexose 
metabolism. Conversely, lipid metabolic process, signal 
transduction, response to stress, and response to 



stimulus were mainly repressed by P. viticola (60, 58, 56 
and 55%, respectively); this included genes encoding 
protein kinases, receptor-like protein kinases, disease re- 
sistance proteins, and PRs. 

The categories of response to stress, secondary metab- 
olism, and DNA metabolism were significantly overre- 
presented in clusters 6 and 7, and the genes in these 
categories were mainly up-regulated following P, viticola 
inoculation (Figure 4d). In particular, expression of pro- 
tein kinase, disease resistance protein, chitinase {Chit), 
stilbene synthase {STS)y resveratrol o-methyltransf erase, 
phenylalanine ammonia-lyase {PAL), and transcription 
factor genes was enhanced after P, viticola inoculation 
in T39-treated plants compared with control plants. 

Clusters of genes modulated by P. viticola exclusively 
in control plants revealed global down-regulation of 
genes related to protein metabolic process (62 and 71% 
in clusters 8 and 9, respectively), response to stress 
(70 and 60%, respectively), and response to stimulus 
(66 and 57%, respectively; Figure 4e). Interestingly, the 
same categories repressed by P. viticola in control 
plants were globally induced in T39-treated plants 
(cluster 10; Figure 4f). In particular, genes related to 
protein metabolic process, response to stress, response 
to stimulus, and transcription were mainly up-regulated 
by P, viticola in T39-treated plants (77, 67, 61, and 
67%, respectively). This active reaction to the pathogen 
in T39-treated plants included expression of protein 
kinases, transcription factors, auxin and JA/ET signals, 
phenylpropanoid biosynthesis, and defence-related 
genes. 
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Figure 4 (See legend on next page.) 
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Validation of RNA-Seq analysis 

To validate the RNA-Seq results, we used real-time RT- 
PCR to analyse the expression levels of 24 differentially 
expressed genes (Figure 5). We selected genes with dif- 
ferent expression profiles and expression levels, 



including genes associated with different functional cat- 
egories (Additional file 11). Although the extent of 
modulation revealed by real-time RT-PCR and RNA-Seq 
may differ [6,48,49], the real-time RT-PCR expression 
profiles in our analysis were in complete agreement with 
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(See figure on previous page.) 

Figure 4 Functional categories of differentially expressed genes grouped into clusters based on their expression profiles. Numbers of 
induced (red) and repressed (green) grapevine genes witliin eacli selected Gene Ontology (GO) functional category are reported for each cluster: 
(A) Genes modulated by Trichodermo horzionum T39 (T39) with comparable expression levels before and after Plosmoporo viticolo inoculation in 
T39-treated plants (cluster 1, solid bars) and genes initially modulated by T39 and further modulated by P. viticolo inoculation in T39-treated 
plants (cluster 2, hatched bars); (B) Genes modulated exclusively by T39 treatment (cluster 3); (C) Genes similarly modulated by P. viticolo in 
control and T39-treated plants with (cluster 4, solid bars) or without (cluster 5, hatched bars) direct modulation by T39 treatment; (D) Genes 
modulated by P. viticolo with greater modulation (> 1.5-fold) in T39-treated than in control plants (ISR-primed genes) with (cluster 6, solid bars) or 
without (cluster 7, hatched bars) direct modulation by T39 treatment; (E) Genes modulated by P. viticolo exclusively in control plants with (cluster 
8, solid bars) or without (cluster 9, hatched bars) direct modulation by T39; (F) Genes modulated by P. viticolo exclusively in T39-treated plants 
(cluster 10). For each cluster of differentially expressed genes, GO functional categories marked by an asterisk were differentially represented to 
the entire grapevine transcriptome, according to GOstat statistical analysis {P< 0.05) [87]. 



the RNA-Seq data for 19 genes. The expression profiles 
generated by real-time RT-PCR and RNA-Seq differed 
for three genes in one treatment and for two genes in all 
treatments. These differences could be due to differences 
in sensitivity, particularly in distinguishing members of 
multigene families. Specificity for RNA-Seq analysis was 
observed in two grapevine chitinases (Table 4) that are 
known to be strongly {Chit la) or weakly {Chit lb) up- 
regulated after P, viticola inoculation [50]. 

We also validated the RNA-Seq method by confirming 
the expression of previously identified markers of plant 
response to Trichoderma spp. and P. viticola (Table 4). 
Our RNA-Seq data confirmed P. v/^/co/<3^- dependent up- 
regulation of PR-l, PR'2, PR'4, and PR-IO genes in sus- 
ceptible grapevines [5,7,11,19,51,52]. Expression profiles 
of PR'2, PR'4, PR'S, PR'IO, and lipoxygenase 9 (IOX-9) 
were consistent with previous real-time RT-PCR data 
obtained during T39-induced resistance [42]. The RNA- 
Seq analysis revealed T39-dependent induction of genes 
that are known to be involved in TrichodermaAnduced 
resistance in other systems [26,29,53,54]. Moreover, 
RNA-Seq analysis confirmed the modulation of genes 
known to be affected by P, viticola in grapevines 
[4-7,9,11,17,52]. 

Discussion 

RNA-Seq analyses of leaf samples collected before and 
24 h after P, viticola inoculation from T39-treated and 
control plants resulted in the identification of 7,024 dif- 
ferentially expressed genes. These genes formed 10 clus- 
ters of different expression profiles, highlighting the 
complex transcriptional reprogramming of T39-induced 
resistance. T39 treatment directly affected the expression 
of grapevine genes and to a greater extent enhanced 
grapevine response to P, viticola inoculation, indicating 
a dual effect of T39. At the sampling time points 
selected in this analysis, we showed that a limited num- 
ber of changes in gene expression were caused by T39 
treatment and that more intense transcriptional repro- 
gramming took place after pathogen inoculation. In par- 
ticular, opposing modulation of genes related to 
response to stress and protein metabolism was observed 



in T39-treated and control plants at 24 h after P, viticola 
inoculation, indicating inhibition of disease-related pro- 
cesses and induction of active defence in T39-treated 
grapevines. Based on the expression profiles, genes dir- 
ectly modulated by T39, as well as genes with reinforced 
or specific modulation in T39-treated plants after patho- 
gen inoculation are strong candidates for activation of 
plant self-protection and consequent inhibition of 
disease-related processes and symptoms development. 

Grapevine processes directly affected by Trichoderma 
harzianum T39 

Analysis of the expression profiles revealed a set of genes 
directly modulated by T39 (182 up- and 161 down-regu- 
lated) and whose expression was not affected by the sub- 
sequent pathogen inoculation (cluster 1). Several genes 
of cluster 1 were involved in signal transduction pro- 
cesses, indicating that they may be related to the initial 
events of recognition of the beneficial microorganism 
and induction of resistance, as recently demonstrated by 
the proteomic analysis of T39-induced resistance [55]. 
Enzymes known to mediate microbial recognition and to 
trigger defence responses in plant species [56] were up- 
regulated by T39, and these included 66 receptor-like 
protein kinase genes, three protein kinase genes, and 
one protein phosphatase gene. In particular, the serine/ 
threonine kinase receptor {STKR) was similar to the 
Arabidopsis FRKl induced by bacterial flagellins [57], 
and a protein kinase was homologous to the Tricho- 
dermaAnduced kinase {TIPK) of cucumber [53]. Further- 
more, the STKR gene (JG391826) is induced by P. 
viticola in resistant grapevines [8], suggesting that it 
plays a crucial role in the activation of specific defence 
processes against downy mildew. 

Our data suggest that resistance induction was 
also mediated by hormone signalling and transcriptional 
reprogramming. Genes related to ET metabolism 
(two 1-aminocyclopropane-l-carboxylate oxidases) were 
induced by T39 treatment, together with two MYB 
(including MYB72), one FAMA, and one NAC transcrip- 
tion factors. Activation of ET metabolism was in agree- 
ment with previous data showing that JA/ET signals are 
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(See figure on previous page.) 

Figure 5 Comparison of RNA-Seq and real-time RT-PCR analyses. Expression profiles of (A) transcription factor FAMA, (B) Trichoderma- 
induced osmotin 2 {0SM2), (C) serine/tlireonine l<inase receptor (STKR), (D) liypotlietical protein {HP), (E) transcription factor of cluster 2 
{WRKY-C2), (F) WRKY transcription factor 1 (WRKYl), (G) chalcone synthase 3 {CHS3), (H) auxin response factor (ARF), (!) constans-like protein (COL), 
(J) mitogen-activated protein kinase [MARK), (K) protein phosphatase 2C {PP2Q, (L) polyphenol oxidase [PRO), (M) chlorophyll a-b binding protein 
(CAB), (N) acidic endochitinase 3 {ChitS), (O) Tr/c/iodermcz-induced osmotin 1 (OSMl), (?) AP2-like ethylene-responsive transcription factor {AP2), 
(Q) invertase (/A/1//), (R) receptor-like serine/threonine-protein kinase {RLK), (S) universal stress A-like protein {UspA), (T) ABC transporter C family 
member {ABC-Q, (U) sucrose synthase (55), (V) ethylene-responsive transcription factor {ERF), (W) LRR receptor-like serine/threonine-protein kinase 
{LK), and (X) WRKY transcription factor of cluster 10 [WRKY-CIO). Blue dotted lines represent expression levels (FPKM) as assessed by RNA-Seq 
analysis and reported as means and standard errors of three biological replicates for each treatment: control (C), Trichodermo horzionum T39- 
treated (T39), Plasmopara wt/co/c/-inoculated control {C+P.v), and P. wt/co/a-inoculated T39-treated {J39+P.v) plants. Histograms represent the 
relative expression levels (fold-change relative to the expression in control plants) as assessed by real-time RT-PCR and reported as means and 
standard errors of three biological replicates for each treatment. 



involved in T39-induced resistance [42,58]. In Arabidop- 
siSy MYB72 plays a crucial role in TrichodermaAnduced 
resistance [26] and is required in early signalling of 
rhizobacteria-induced resistance in Arabidopsis [59]. 
The NAC gene is involved during oxidative stress [60], 
and response to oxidative stress after T39 treatment was 
suggested by modulation of two peroxidases, three gluta- 
thione S -transferases {GSTs), and one thioredoxin. More- 
over, defence-related genes (two STSs, a Chit, and a 
cellulose synthase) were also directly induced, indicating 
that the maximum expression level of these genes was 
probably reached after T39 treatment and that this was 
sufficient to contribute to defence against subsequent 
P. viticola inoculation. Other defence-related genes were 
pre-induced by T39 and further induced after P, viticola 
inoculation (cluster 2), indicating reinforcement of 
grapevine defence after pathogen inoculation. 

Enhanced reaction of T39-treated plants to 
Plasmopara viticola 

Plant resistance induced by beneficial microorganisms 
has been associated with faster and/or stronger activa- 
tion of defence responses after pathogen inoculation 
[37]. P. v/^/oc/<3^-responsive genes with enhanced expres- 
sion in T39-treated plants were clustered to distinguish 
between those directly modulated (cluster 6) and those 
not modulated (cluster 7) by T39 treatment. These expres- 
sion profiles provide strong support for the view that Tri- 
choderma spp. may have a dual effect: it directly 
stimulates induction of some genes and ftirther reinforces 
modulation of these and other genes after pathogen inocu- 
lation [42] . The dual effect was also reported for defence- 
gene modulation during resistance induced by sulfated 
laminarin [17] and for phytoalexin accumulation during 
resistance induced by Rheum palmatum extracts [24]. 

Enhanced reaction of T39-treated plants included in- 
duction of genes related to response to stimulus and re- 
sponse to stress categories, suggesting improvement of 
signalling pathways and activation of defence reactions 
in response to pathogen inoculation. In particular, 50 
receptor-like kinase genes (35 in cluster 6 and 15 in 



cluster 7) and six protein kinase genes (one in cluster 6 
and five in cluster 7) of signal transduction were primed 
in T39-treated plants. In agreement with our results, 
priming was associated with increased expression of 
mitogen-activated protein kinases {MPK3 and MPK6) in 
Arabidopsis plants [43]. Up-regulation of MPK genes 
was also associated with the reaction of resistant grape- 
vines to P, viticola [4,6], and post-translational modifica- 
tion may be the additional mechanism for activating 
MPKs in response to pathogens [20]. Moreover, expres- 
sion of stress-related genes was enhanced in inoculated 
T39-treated plants, and these genes included PR'2, PR'4, 
PR'S, PR'IO, Chits, and osmotin [OSMl). Similar pro- 
files were observed for Arabidopsis PR genes during T. 
asperelloides T203 and Pseudomonas syringae interac- 
tions [30], indicating common mechanisms of Tricho- 
dermaAnduced resistance against biotrophs. Primed 
profiles were found in resistance-related genes (12 genes 
in cluster 6 and eight in cluster 7). Peronospora parasi- 
tica resistance genes {RPP) must be expressed at optimal 
levels to function against downy mildew in Arabidopsis 
[61], and specific profiles of grapevine RPP genes suggest 
fine tuning during T39-induced resistance. 

A key role of transcriptional regulation and secondary 
metabolic processes in T39-induced resistance was also 
indicated by the priming profiles of four MYB, two 
WRKY, and 12 ET- responsive transcription factors (e.g., 
AP2), and of PAL and STSs genes. Moreover, the cat- 
egory DNA metabolic process was significantly overre- 
presented in ISR-primed genes. In particular, expression 
of the histone-lysine N-methyltransferases SUVR4 and 
ATX2 was enhanced after P. viticola inoculation, reflect- 
ing possible involvement of epigenetic modifications 
[62-64] in T39-induced resistance. 

Genes with enhanced expression in T39-treated plants 
have been previously related to defence against downy 
mildew in resistant grapevines. This was the case for PR 
genes {PR-2, PR-4, PR- 10, and OSMl), invertase genes 
(INVl and INV2), and genes related to secondary meta- 
bolism {PAL, STSs, STKR, copper-containing amine 
oxidase, polyphenol oxidase, and three resveratrol 
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Table 4 Expression levels of grapevine genes known to be modulated during Trichoderma spp.-induced resistance and 
in response to Plasmopara viticola inoculation 
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51.3 


Repressed by P. viticola [9] 


B-amylase {AMY) 


fgenesh.W78X223906.2_l 


3.7 


3.2 


1.4 


1.2 


Repressed by P. viticola [1 0] 


B-amylase {AMY) 


glimmer.W78X087147.9_2 


57.8 


52.7 


11.8 


12.0 


Repressed by P. viticola [1 0] 


sucrose pliospliate syntliase {SPS) 


fgenesh.W78X130197.5_l 


111.8 


78.1 


32.2 


28.8 


Repressed by P. viticola [1 0] 


Pliosplioglucan, water dil<inase {PWD) 


fgenesh.W78X092604.4_l 


23.4 


17.0 


12.1 


10.3 


Repressed by P. viticola [1 0] 


Glucose-1 -pliospliate adenylyltransferase 
{GPAD 


sim4.W78Xl 53373.1 1_4 


103.6 


65.5 


39.2 


29.3 


Repressed by P. viticola [1 0] 


Tlirealase 


fgenesh.W78X152339.9_5 


1.0 


1.3 


2.7 


2.8 


Induced by P. viticola [10] 


a,a-trelialose-pliospliate syntliase {TPS) 


glimmer.W78X200879.5_4 


0.3 


0.2 


0.1 


0.1 


Repressed by P. viticola [1 0] 


Plienylalanine ammonia lyase {PAL) 


glimmer.W78X178553.4_4 


0.0 


0.2 


0.2 


0.3 


Induced by P. viticola [17,50,94] 
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Table 4 Expression levels of grapevine genes known to be modulated during Trichoderma spp.-induced resistance and 
in response to Plasmopara viticola inoculation (Continued) 
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n 1 

U. 1 
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1 4 


2 2 
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1 0 
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InHi ir"pH \~wi P \/ifimiln TAI 
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(3-glucosidase {BG) 


fgenesh.W78X148924.9_l 


0.1 


0.2 


2.0 


0.7 


Induced by P. viticolo [6] 


TMV response-related protein [TMVR) 


fgenesh.W78Xl 06668.41 _2 


13.5 


7.1 


8.6 


6.7 


Induced by P. viticolo [6] 


Polygalacturonase-inhibiting protein {PGIPf 


glimmer.W78X089275.7_l 










Induced by P. viticolo [6] 


WRKY transcription factor 21 {WRKY21) 


fgenesh.W78X121610.8_2 


32.4 


19.9 


6.5 


6.9 


Induced by P. viticolo [6] 


Etiiylene-regulated transcript 2 {ERT2f 


glimmer.W78X175653.7_2 










Induced by P. viticolo [6] 


Rapidly elicited Avr9/Cf-9 proteins {Avr/Cf9) 


fgenesh.W78X009646.1 2_6 


71.5 


70.8 


17.9 


21.7 


Modulated by P. viticolo [5] 


Rapidly elicited Avr9/Cf-9 proteins {Avr/Cf9) 


fgenesh.W78X263622.1 9_2 


236.7 


117.8 


27.1 


22.3 


Induced by P. viticolo [4] 


Harpin-induced protein 1 (Hinl) 


g 1 i mmer.W78X225495.6_2 


53.7 


63.2 


12.1 


16.0 


Induced by P. viticolo [4] 


Hypersensitivity-related 203 J {Hsr203j) 


fgenesh.W78X137613.4_l 


52.3 


45.0 


7.7 


9.1 


Induced by P. viticolo [4,12] 



^ Grapevine genes of the Pinot Noir grapevine genome [77] Release 3 [78]. 

^ Gene expression levels calculated as fragments per kilobase of transcript per million fragments mapped (FPKM) In control (C), Trichoderma harzianum 
T39-treated (T39), Plasmopara v/f/co/a-lnoculated control {C+P.v.), and P. wf/co/a-lnoculated T39-treated (T39+P.V.) plants. 

References of previous analyses. 
^ Genes not Included In the set of differentially expressed genes. 



O-methyltransferases), for which up-regulation was 
greater in resistant than in susceptible grapevines upon 
P. viticola inoculation [4,8,11-14]. The expression pro- 
files of these marker genes suggest that the defence pro- 
cesses usually activated against downy mildew in 
resistant grapevines are partially stimulated in susceptible 
plants by T39-induced resistance. 

Disease-related processes inhibited in T39-treated plants 

Another important aspect of T39-induced resistance 
was evidenced by genes modulated by P, viticola exclu- 
sively in control plants (clusters 8 and 9). These genes 
were mainly down-regulated (57 and 55%, respectively), 
and they reflect exploitation of cellular resources and/ 
or suppression of defence responses during the com- 
patible interaction. Down-regulation caused by P. viti- 
cola in control plants involved categories of response 
to stress and primary metabolic processes. Many signal- 
ling components (kinase, phosphatase, calmodulin, and 
calcium signalling), transcription factors {WRKY and 
MYB), and disease resistance proteins were repressed, 
supporting the view that suppression of endogenous 
signals is required to establish the compatible inter- 
action [5]. Three ABC transporters (e.g., ABC-C) were 
also repressed, and suppression of some ABC transpor- 
ters increases the susceptibility to oomycete pathogens 
[65]. Specific alteration of carbohydrate metabolism by 



P, viticola in control plants was highlighted by modula- 
tion of glucosidase, galactosidase, mannosidase, and su- 
crose synthase {SS) genes. Moreover, repressed genes of 
cluster 9 were classified into the categories of energy 
metabolism (a phytochrome C, two malic enzymes, and 
a ribokinase) and defence response (five chitinases, 
three glucanases, three superoxide dismutases, and a 
callose synthase), reflecting disease-related process 
employed by P, viticola only in control plants. P, viti- 
cola might need to actively suppress plant defences 
during leaf colonisation through microbial effectors, as 
demonstrated in other oomycetes [66,67]. Particularly, 
hemibiotrophic and biotrophic species establish intim- 
ate associations with plants [66,67]. To establish infec- 
tion, these pathogens must suppress the plant defence 
and manipulate the host metabolism by microbial effec- 
tors (virulence factors) that are translocated inside the 
plant cell or secreted into the extracellular space within 
plant tissue [66,68]. P, viticola effector genes have been 
recently described [69], and modulation of host cell 
defences through virulence factors in susceptible grape- 
vines has been indicated by histochemical [70] and 
transcriptomic [4] analyses. In agreement with the 
phenotypic observations, the specific modulation of 
grapevine genes in inoculated control but not in inocu- 
lated T39-treated plants indicates that T39-induced 
resistance acts by inhibiting some disease-related 
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processes and probably by interfering with some pathogen- 
induced processes. 

Specific transcriptional response of T39-treated plants to 
Plasmopam viticola 

In contrast to broad down-regulation of genes in control 
plants (clusters 8 and 9), genes specifically modulated by 
P. viticola in T39-treated plants (cluster 10) were mainly 
induced (63%). These opposing reactions to P. viticola 
are particularly evident in genes related to protein me- 
tabolism, response to stimulus, and response to stress, 
which were mainly induced in T39-treated plants and 
mainly repressed in control plants. Up-regulation of 
genes associated with these categories has been observed 
in resistant grapevines [4,13,71], indicating that the de- 
fence processes of resistant genotypes could be partially 
activated in susceptible varieties by T39-induced resist- 
ance. In particular, all NBS -encoding resistance (NBS-R) 
genes modulated by P, viticola in control plants (cluster 
9) were repressed, whereas those with ISR-responsive 
specific profiles (cluster 10) or ISR-primed profiles (clus- 
ter 7) were mainly up-regulated. Opposing modulation 
of NBS'R genes probably reflects suppression of plant 
defence in control plants and activation of defence 
responses in T39-treated plants. Interestingly, the NBS- 
R genes of clusters 7 and 10 belong mainly to the Va 
component genome of grapevine [71], indicating 
subgenome-dependent regulation of gene expression 
[72] in grapevine. 

Defence signals specifically activated in T39-treated 
plants included those mediated by auxin (two auxin 
transporters, two auxin-induced proteins, and two 
indole-3-acetic acid amido synthetases), ET (ACC oxi- 
dase and five ERF transcription factors), and J A (three 
lipoxygenases and two fatty acid desaturases). The role 
of JA/ET signalling pathways in T39-induced resistance 
has also been demonstrated by phytohormone-affected 
Arabidopsis mutants [58] and by expression analysis of 
grapevine marker genes [42]. ISR is commonly regulated 
by JA/ET-dependent signalling pathways, and it is espe- 
cially active against pathogens deterred by defences that 
are controlled by JA and ET [33,36]. The auxin response 
pathway is connected to the SA and JA/ET signalling 
networks [33], and crosstalk between hormonal path- 
ways [36] enables the fine tuning of defence mechanisms 
so that the plant can tailor its response to the specific 
invader [33]. ET exerts its resistance-stimulating activity 
in concert with JA [33], and J A pathways are involved in 
the reaction to P, viticola in resistant grapevines [4,6]. 
Thus, enhancement of JA/ET signals in T39-treated 
plants supports the view that increased resistance to 
downy mildew is mediated by partial activation of extant 
defence mechanisms normally activated in resistant gen- 
otypes. Reaction to the pathogen was also mediated by 



specific up-regulation of 59 receptor kinases, 10 protein 
kinases, two bHLHs genes, one MYB gene, and the 
NPRl.l gene. 

Our results also suggest that the cell redox balance is 
altered in T39-treated plants after pathogen inoculation. 
The reaction of T39-treated plants to pathogen inocula- 
tion included the induction four peroxidases and a GST, 
Antioxidant enzymes are often induced in response to 
pathogens, and alteration of oxidative-stress metabolism 
has a prominent role in the T39-induced resistance of 
grapevine to downy mildew [55]. Peroxidases play sev- 
eral important roles in pathogen resistance by contribut- 
ing to the production of reactive oxygen species, the 
reinforcement of cell walls, and the production of phy- 
toalexins. Accumulation of stilbene phytoalexin is one of 
the most important defence processes activated by re- 
sistant grapevines in response to P, viticola [1]; genes of 
phenylpropanoid biosynthesis (flavanone 3-dioxygenase, 
laccase, and dihydroflavonol-4-reductase) were specific- 
ally induced in T39-treated plants, confirming activation 
of pathways known in resistant genotypes. However, 
additional defence mechanisms against P. viticola are 
activated in resistant genotypes. In particular, HR-related 
genes {Avr/Cf9, Hinl, and Hsr203j) were not induced 
and localised HR necrosis was not observed in T39- 
treated plants. 

Common transcriptional response of control and 
T39-treated plants to Plasmopara viticola 

Although specific transcriptional reprogramming of 
T39-treated plants was observed, 3,454 genes had com- 
parable expression levels in control and T39-treated 
plants after P, viticola inoculation (clusters 4 and 5). The 
pathogen-responsive processes not affected by resistance 
induction were mainly related to primary metabolism 
and signal transduction. In particular, expression profiles 
of genes related to starch metabolism (up-regulation of 
a-amylase and sugar transporters, down-regulation of 
|3-amylase, glucose- 1 -phosphate adenylyltransferases, 
sucrose phosphate synthase, and phosphoglucan water 
dikinase) indicated the source-to-sink transition of P, 
viticola-infected leaves [10]. Likewise, genes related to 
photosynthesis (two quinone oxidoreductases, a chlor- 
ophyllide oxygenase, a protochlorophyllide transporter, 
a chlorophyllase-2 and chlorophyll a-b binding pro- 
teins) and to the Calvin cycle (two rubisco genes and 
a phosphoglycerate kinase gene) were similarly modu- 
lated in inoculated control and T39-treated plants, 
possibly reflecting the establishment of a compatible 
interaction. P, viticola inoculation also resulted in the 
down-regulation of genes involved in the signal trans- 
duction processes (20 receptor-like protein kinase 
genes, 11 protein kinase genes, and five protein phos- 
phatase genes) and defence response (23 probable 
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disease resistance genes, six unspecified PR genes, a 
GST gene and a thaumatin-like gene [73]), indicating 
a pathogen-dependent suppression of the host reac- 
tion mechanisms. 

Other markers of P, viticola infection in susceptible 
grapevines showed comparable modulation in T39-treated 
and control plants, such as enzymatic resistance protein 
[5], zeaxanthin epoxidase [9], catalase [52], isoforms of 
chalcone synthase [7], and dihydroflavonol 4-reductase 
[11]. These results suggested that transcriptional changes 
associated with the compatible interaction [5,7,10,13] were 
not completely inhibited in T39-treated plants, which is 
consistent with the observation that downy mildew symp- 
toms were reduced but not completely blocked by T39 
treatment. 

Conclusions 

The transcriptome analysis reported here represents a 
major contribution to the characterization of induced re- 
sistance mechanisms in a non-model plant. We used the 
RNA-Seq approach to characterize the transcriptional 
changes associated with the reduction of downy mildew 
symptoms in T39-treated grapevines and thereby estab- 
lished a foundation for a more detailed time-course ana- 
lysis of induced resistance in grapevines. Our data show 
that T39 directly activates microbial recognition 
mechanisms in absence of pathogen infection and 
enhances the expression of defence processes after P, 
viticola inoculation. Reduction of downy mildew symp- 
toms is related to inhibition of disease-related processes 
and activation of defence mechanisms after P. viticola in- 
oculation of T39-treated plants. T39-induced resistance 
is associated with enhanced expression of specific genes 
related to resistance against downy mildew in wild grape- 
vines, indicating that induced resistance can partially 
mimic defence processes of resistant genotypes. In 
addition, the genes identified in this work represent an 
important source of markers of resistance induction and 
can be used to select beneficial microorganisms with an 
improved ability to induce plant resistance. These mar- 
kers can also be used to clarify how environmental con- 
ditions affect induced resistance. 

Methods 

Biological materials 

Two-year-old plants of the susceptible grapevine {V, 
vinifera) cultivar Pinot Noir grafted onto Kober 5BB 
were individually planted in 2.5 L pots containing a mix- 
ture of peat and pumice (3:1). Plants were grown for two 
months in a greenhouse at 25 ± 1°C with a photoperiod 
of 16 h light and relative humidity (RH) of 70 ± 10%. T. 
harzianum T39 (Trichodex, Makhteshim Ltd., Israel) 
was applied to grapevine leaves at 8 g L"^ in water, corre- 
sponding to a conidial suspension of 10^ colony-forming 



units (cfu) mL'^. A P. viticola isolate was collected from 
an untreated vineyard in the Trentino region (northern 
Italy) and maintained by subsequent inoculations on 
V, vinifera Pinot Noir plants under controlled green- 
house conditions. Plants with oil spot symptoms were 
incubated overnight in the dark at 99-100% RH. 
P. viticola sporangia were collected by washing the ab- 
axial surfaces bearing freshly sporulating lesions with 
cold distilled water, and the inoculum concentration was 
adjusted to 10^ sporangia mL'^ with a haemocytometer. 

Bioassay of induced resistance in grapevine 

All leaves of each shoot were treated three times with 
T39 (at one, two and three days before pathogen inocu- 
lation) to induce the greatest phenotypic resistance re- 
sponse [21]. Leaves directly treated with the resistance 
inducer were analysed to study local effects of T39- 
induced resistance, and untreated plants were used as 
control. The suspensions were applied to the abaxial and 
adaxial leaf surfaces using a compressed air hand sprayer 
(20-30 mL for each plant). 

One day after the last T39 treatment, the abaxial sur- 
faces of all leaves of each plant were inoculated with a 
P, viticola suspension (10^ sporangia mL'^) using a 
compressed-air hand sprayer (20-30 mL for each 
plant). All plants were then incubated overnight in the 
dark at 25 ± 1°C with 99-100% RH and then kept 
under controlled greenhouse conditions. 

Six days after inoculation with the pathogen, all plants 
were incubated overnight in the dark at 25°C with 99- 
100% RH. Disease severity was visually assessed as the 
percentage of abaxial leaf area covered by sporulation 
[74]. Each treatment was carried out on six replicates 
(plants). An ANOVA analysis was performed using the 
Statistica 9 software (StatSoft, Tulsa, OK) followed by 
Tukey s test to detect significant differences {P < 0.05) in 
disease severity between treatments. 

Sample collection and RNA isolation 

Leaf samples were collected from control and T39- 
treated plants immediately before and 24 h after P, viti- 
cola inoculation. This time point was chosen because it 
is associated with leaf colonization by primary hyphae 
[24,55,75] and with modulation of defence-related genes 
[4,17,19,42] for the establishment of defence responses 
in resistant genotypes [4] . Four treatments were analysed 
by RNA-Seq: control (C), T39-treated (T39), P. viticola- 
inoculated control (C+P.v.), and P, viticola-inoculsited 
T39-treated (T39+P.V.) plants. For each treatment, leaf 
samples from three replicates (plants) were collected at 
each time point, i.e., three plants per treatment were 
sampled before inoculation and three plants per treat- 
ment were sampled 24 h after inoculation; different 
plants were sampled at each time point to avoid the 
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effects of wounding stress. Each sample comprised three 
leaves taken from the same plant, and only leaves of the 
4^^-5^^ node from the top of the shoot were collected to 
avoid ontogenic resistance effects [42]. Samples were im- 
mediately frozen in liquid N2 and stored at -80°C. Total 
RNA was extracted using the Spectrum Plant total RNA 
kit ( Sigma- Aldrich, St. Louis, MO) and quantified using 
the Nanodrop 8000 (Thermo Fisher Scientific, Wilmington, 
DE). The quality of the RNA extracts was checked 
using the Bioanalyzer 2100 (Agilent Technologies, 
Santa Clara, CA). 

Library construction and lllumina sequencing 

Three biological replicates for each treatment were sub- 
jected to RNA-Seq library construction, and each library 
was sequenced twice in separate lanes (sequencing repli- 
cate). In a preliminary experiment, a technical replicate 
of one sample was run, and the results showed that the 
library preparation had a high level of technical reprodu- 
cibility (R^ = 0.93, Pearson s correlation = 0.96). Libraries 
were prepared using the TruSeq SBS v5 protocol 
(lllumina, San Diego, CA), and paired-end reads of 100 
nucleotides were obtained using an lllumina HiSeq 2000 
at Fasteris (Plan-les-Ouates, Switzerland). Briefly, poly 
(A) -containing mRNA was isolated from 10 (ig of total 
RNA in two rounds of purification using poly-T oligo- 
attached magnetic beads. Purified mRNA was fragmen- 
ted using Zn-catalysed hydrolysis and converted into 
double-stranded cDNA by random priming. Following 
end repair, indexed adapters were ligated and cDNA 
fragments of 200 ± 25 bp were purified. Purified cDNA 
was amplified by PCR and validated by Sanger sequen- 
cing, after which mRNA-Seq libraries were multiplexed 
(six libraries per lane) and sequenced according to the 
manufacturers instructions. The sequences have been 
deposited at the Sequence Read Archive of the National 
Center for Biotechnology [76] under BioProject number 
PRJNA168987. 

Mapping of sequenced reads and assessment of gene 
expression 

For sequence quality control, in-house python scripts 
were used for quality trimming and sequence filtering: 
sequencing adapters, k-mers, and bases with a Phred 
quality score lower than 30 were removed from the read 
ends, and reads shorter than 50 nucleotides were then 
discarded. Filtered reads were mapped to the Pinot Noir 
grapevine genome [77] Release 3 [78] using TopHat 
1.2.0 release with default settings [79]; up to 100 hits for 
each read were allowed to improve gene expression 
evaluation of members of multigene families [45]. Gene 
expression values were determined using Cufflinks 1.0.3 
release [46], and the FPKM values were calculated for 
each transcript. Cufflinks default settings were adopted 



and gene FPKM values were computed by summing the 
FPKM values of different transcripts of the same gene 
[46]. The option for fragment bias correction was ap- 
plied [80], and the RABT assembly option was used to 
improve the identification of novel genes and transcripts 
[81]. Sequencing replicates were first analysed separately, 
and then the read counts of two sequencing replicates 
were summed and FPKM values recalculated for each 
biological replicate. 

Differential gene expression analysis 

The number of reads falling into each transcript was 
estimated from the FPKM values according to the for- 
mula reported by Zenoni et al. [49], and gene counts 
were computed by summing counts of different tran- 
scripts of the same gene [46]. Differentially expressed 
genes were identified by the DESeq package [47], with 
an FDR of the Benjamini-Hochberg multiple tests of 5% 
(P<0.05) based on read counts and a minimum fold- 
change of two in at least one pairwise comparison: T39 
versus C, C+P.v. versus C, or T39+P.V. versus C. Fold- 
changes were calculated from FPKM expression values 
with the minimum expression value (10'^ FPKM) 
imposed on non-expressed genes. Genes modulated by 
pathogen inoculation exclusively in T39-treated plants 
were classified as ISR-responsive specific genes [38]. 
Among genes modulated by P. viticola in control and 
T39-treated plants, genes that showed an up- or down- 
regulation greater than 1.5 -fold in inoculated T39- 
treated versus inoculated control plants were selected as 
ISR-primed genes. This priming cut-off value is based 
on previous expression data of Arabidopsis genes that 
showed a primed expression profile after pathogen at- 
tack [38]. ISR-primed genes were further distinguished 
as genes directly modulated (cluster 6) or not (cluster 7) 
by T39 treatment. Principal component analysis (PCA) 
was performed with Statistica 9 software (StatSoft, Tulsa, 
OK) on FPKM expression values of all grapevine genes, 
and the expression profiles of differentially expressed 
genes were determined by cluster analysis using the 
T-MeV 4.8.1 software [82]. 

Gene annotation 

The grapevine genes were automatically annotated using 
the ARGOT2 function prediction tool [83,84], which 
was able to annotate grapevine protein sequences by 
blastp search against the UniProtKB database [85] 
(downloaded from January 2012) and by association of 
Gene Ontology biological process terms [86] weighted 
according to both semantic similarity relations and asso- 
ciated scores [83,84]. Annotation of differentially 
expressed genes and their distribution over 15 selected 
functional categories was compared with the entire 
grapevine transcriptome using GOstat statistical analysis 
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[87]. For manual annotation, grapevine proteins were 
aligned against the Swiss Prot [88] database (downloaded 
from February 2012) using blastp, and the best hits 
(E-value lower than lE'^) were used to select protein 
descriptions. The corresponding genes of the grapevine 
PN40024 genotype [89] were identified by blastp search 
against the 12 x release [90] with an imposed E- value 
lower than lE"^^. 

Gene expression analysis by quantitative real-time RT-PCR 

Total RNA was treated with DNase I (Invitrogen, 
Carlsbad, CA), and the first-strand cDNA was synthe- 
sized from 1.0 (ig of total RNA using Superscript III 
(Invitrogen) and oligo-dT. Reactions were carried out 
with Platinum SYBR Green qPCR SuperMix-UDG (Invi- 
trogen) and specific primers (Additional file 11) using 
the Light Cycler 480 (Roche Diagnostics, Germany). The 
PCR conditions were: 50°C for 2 min and 95°C for 2 min 
as initial steps, followed by 50 cycles of 95°C for 15 s 
and 60°C for 1 min. Each sample was examined in three 
technical replicates, and dissociation curves were ana- 
lysed to verify the specificity of each amplification reac- 
tion. Light Cycler 480 SVl.5.0 software (Roche) was 
used to extract Ct values using the second derivative cal- 
culation [91], and LinReg software was used to calculate 
reaction efficiencies [92]. The relative expression of each 
gene was then calculated according to the Pfaffl equation 
[93] using control plants as the calibrator. The actin 
gene (glimmer. VV78X1 149 14.6_2) was used as the con- 
stitutive gene for normalization [4,52] because its ex- 
pression was not significantly affected by the treatments 
(Additional file 12), and comparable results were 
obtained with VATP16 [94]. Mean expression and stand- 
ard error of three biological replicates were calculated 
for each treatment. 

Additional files 



Additional file 1: RNA-Seq sequencing and read mapping of each 
sequencing replicate. Numbers of RNA-Seq reads passir^g the quality 
check and mapping to the grapevine genome are reported for each 
sequencing replicate (named A and B) of each biological replicate 
(numbered from 1 to 3) for control (C), Trichoderma harzianum T39- 
treated (T39), Plasmopara wr/co/a-inoculated control (C+P.v.), and P. 
wT/co/a-inoculated T39-treated (J39+P.V.) plants. 

Additional file 2: RNA-Seq sequencing and coverage of the 
grapevine transcriptome for each biological and sequencing 
replicate. Total bases (Mbp) sequenced by RNA-Seq analysis and 
coverage of the grapevine transcriptome are reported for each 
sequencing replicate (named A and B) of each biological replicate 
(numbered from 1 to 3) for control (C), Trichoderma harzianum T39- 
treated (T39), Plasmopara wT/co/a-inoculated control {C+P.v.), and P 
Wr/co/c/-inoculated T39-treated (J39+P.V.) plants. 

Additional file 3: Distribution of RNA-Seq sequences within the 
grapevine genome. (A) Proportions of reads mapping to unique 
locations (unique reads, yellow), reads mapping to multiple locations 
with 2-100 matches (multi-reads, red), and reads not mapping or 



ambiguously mapping to more than 100 locations (unused reads, white) 
to the grapevine genome. The mean number of sequenced reads and 
standard errors for three biological replicates are presented for each 
treatment. Percentages (%) are calculated with respect to the total 
filtered reads. Unique reads and multi-reads were used for evaluating 
gene expression. (B) Proportions of reads mapping to grapevine genes 
(green) and to intergenic regions (red) of the grapevine genome. The 
mean number of sequenced reads and standard errors for three 
biological replicates are presented for each treatment. Percentages (%) 
are calculated with respect to the total mapping reads. Grapevine 
treatments: control (C), Trichoderma harzianum T39-treated (T39), 
Plasmopara wT/co/a-inoculated control {C+P.v.), and P wT/co/a-inoculated 
T39-treated (139+P.v) plants. 

Additional file 4: RNA-Seq reads mapping to grapevine genes. 

Numbers of RNA-Seq reads mapping to grapevine genes are reported for 
each sequencing replicate (named A and B) of each biological replicate 
(numbered from 1 to 3) for control (C), Trichoderma harzianum T39- 
treated (T39), Plasmopara w'r/co/a-inoculated control {C+P.v), and P. 
w'r/co/a-inoculated T39-treated (J39+P.V.) plants. 

Additional file 5: Expressed grapevine genes and novel genes 
identified by RNA-Seq analysis. Numbers of expressed grapevine 
genes and novel genes with respect to the grapevine annotation are 
reported for each sequencing replicate (named A and B) of each 
biological replicate (numbered from 1 to 3) for control (C), Trichoderma 
harzianum T39-treated (T39), Plasmopara w'r/'co/a-inoculated control (C+P. 
v.), and P. wT/'co/a-inoculated T39-treated (T39+P.^.) plants 

Additional file 6: Expression levels of grapevine genes. Gene 
expression values (FPKM) and standard errors are reported for all (33,514) 
Pinot Noir grapevine predicated genes [77] Release 3 [78] in control (C), 
Trichoderma harzianum T39-treated (T39), Plasmopara w'r/'co/a-inoculated 
control {C+P.v), and P. wT/co/a-inoculated T39-treated (T39+P.^.) plants. 
Expression values in each biological replicate (numbered from 1 to 3) 
and each sequencing replicate (named A and B) are also listed. 

Additional file 7: Correlations between sequencing replicates of 
RNA-Seq analysis. Comparison of the expression levels of all grapevine 
genes, expressed as fragments per kilobase of transcript per million 
fragments mapped (FPKM), in the two sequencing replicates (named A 
and B) of (A) control (C) biological replicate no. 1; (B) C biological 
replicate no. 2 (B); (C) C biological replicate no. 3; (D) Trichoderma 
harzianum T39-treated (T39) biological replicate no. 1; (E) T39-treated 
biological replicate no. 2; (F) T39-treated biological replicate no. 3; (G) 
Plasmopara w'r/co/a-inoculated control {C+P.v.) biological replicate no. 1; 
(H) C+P.v. biological replicate no. 2; (I) C+P.v. biological replicate no. 3; (J) 
P wT/co/a-inoculated T39-treated plants (T39+P.v) biological replicate no. 
1; (K) T39+P.V biological replicate no. 2; and (L) T39+P.V biological 
replicate no. 3. 

Additional file 8: Principal component analysis of grapevine 
treatments. Principal component analysis (PCA) is based on the 
expression values (FPKM) of all grapevine genes for each sequencing 
replicate (named A and B) of each biological replicate (numbered from 1 
to 3) for control (C), Trichoderma harzianum T39-treated (T39), Plasmopara 
w'r/co/a-inoculated control {C+P.v), and P. w'r/'co/c?-inoculated T39-treated 
(Y39+P.V.) plants. 

Additional file 9: Pearson's correlation coefficients between 
sequencing and biological replicates of RNA-Seq analysis. Pearson's 
correlation coefficients are based on the gene expression values (FPKM) 
of all grapevine genes in two sequencing replicates (A and B) of each 
biological replicate (numbered from 1 to 3) for control (C), Trichoderma 
harzianum T39-treated (T39), Plasmopara wt/co/o-inoculated control (C+P. 
v.), and P. wT/'co/a-inoculated T39-treated (T39+P.V.) plants. 

Additional file 10: Expression levels, clustering, and annotation 
results of differentially expressed genes. Gene expression values 
(FPKM) and standard errors are reported for differentially expressed 
(7,024) genes in control (C), Trichoderma harzianum T39-treated (T39), 
Plasmopara w'r/co/a-inoculated control {C+P.v), and P w'r/'co/ci-inoculated 
T39-treated (J39+P.V.) plants. Differentially expressed genes were 
identified by the DESeq package with a false discovery rate (FDR) of 5% 
and a fold-change greater than two in at least one pairwise comparison. 
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Fold-changes, clustering results, and functional annotation are reported 
for each gene. 

Additional file 11: Primer sequences of the grapevine genes 
analysed by real-time RT-PCR. Forward and reverse primer sequences 
of real-time RT-PCR analysis are reported for 24 selected grapevine genes. 
Grapevine /\cr/n (TC81781) [4, 52] and VATP16 (XM_002269086.1) [94] 
were used as constitutive genes for normalising the real-time RT-PCR 
data. 

Additional file 12: Expression levels of grapevine genes that are 
members of the actin gene family. Gene expression values (FPKM) and 
standard errors are reported for grapevine genes belonging to the Actin 
gene family in control (C), Trichoderma harzianum T39-treated (T39), 
Plasmopara w'r/co/ci-inoculated control (C+P.v.), and P. w'r/co/a-inoculated 
T39-treated (T39+P.v) plants. 
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